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The Electrostatic Ion Beam Trap : a mass spectrometer of infinite mass range 

Alexandre VallettcQ and P. Indclicatcjll 

Laboratoire Kastler Brossel, Ecole Normale Superieure, 
CNRS, Universite Pierre et Marie Curie ~ Paris 6, 
Case 74 ; 4, place Jussieu, 75252 Paris CEDEX 05, France 
(Dated: April 22nd, 2011) 

We study the ions dynamics inside an Electrostatic Ion Beam Trap (EIBT) and show that the 
stability of the trapping is ruled by a Hill's equation. This unexpectedly demonstrates that an EIBT 
works like a quadrupole trap. The parallelism between these two kinds of traps is illustrated by 
comparing experimental and theoretical stability diagrams of the EIBT. The main difference with 
quadrupole traps is that the stability depends only on the ratio of the acceleration and trapping 
electrostatic potentials, not on the mass nor the charge of the ions. Our model is confirmed by 
the experimental observation of parametric resonances in the EIBT, where ions excited to their 
oscillation frequency are ejected out of the trap. These frequencies are proportional to the square 
root of the charge/mass ratio of each trapped species, while the stability is independent of the mass 
and charge. The EIBT can thus be used as a mass spectrometer of infinite mass range. 
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Electrostatic Ion Beam Trap (EIBT) [ij are taking 
an important place in between very low-energy charged- 
particles storage devices, such as quadrupole and Pen- 
ning traps [2I and high energy storage rings Q • With the 
ConeTrap [J|, electrostatic rings [5| and the Mini-Ring 
, they form a new family of trap operating at energies 
of a few keV. They are used for atomic and molecular 
metastable-states studies, molecular fragmentation and 
photodissociation (see, e.g., for a review). Beyond pro- 
viding trapping of energetic particles in a well defined di- 
rection, these traps have many interesting features : they 
are small, relatively inexpensive, easy to setup and ope- 
rate and have a field-free region where ions move freely 
and where measurements can easily be performed. They 
can even be used as Time Of Flight (TOE) mass spec- 
trometers [1| or cooled at cryogenic temperatures @. 

Despite these interesting features, all the published 
theoretical models describing EIBT are based on one di- 
mensional approximations, neglecting the radial motion. 
This leads to very inaccurate predictions of the trap sta- 
bility and operating domain, and restrict severely their 
flexibility, as finding reasonable working points requires 
lengthy and tedious experimental exploration. This may 
partly be explained by the lack of an analytical formula 
for the electrostatic potential inside these traps leading 
to the dilemma of choosing between a simplistic analyti- 
cal model and a heavy numerical treatment unsuited to 
explore the huge space of parameters. Usual beam simu- 
lation codes fail to produce good results as the numerical 
inaccuracies at the ion turning points lead to energy non 
conservation reaching a few lOOeV over a few tens of os- 
cillations, which render the simulations useless. Here we 
solve completely the problem, using methods developed 
for radio-frequency quadrupole traps. We show that the 
radial dynamic is ruled by a Hill's equation, a particular 
case of the Mathieu equation that describes quadrupole 
traps. This model yields very accurate predictions of the 
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Figure 1. (Color online) Overview of the experimental setup. 
Five potentials are applied to the electrodes (represented by 
the striped rectangles) , the other are grounded. The injection 
of the ion bunch (orange) is performed when all the electrodes 
on one side of the trap are grounded. The potentials are raised 
before the bunch has time to come back. The faraday cup 
(green) is linked to an oscilloscope via a charge amplifier. 
A sinusoidal voltage can be applied to an initially grounded 
electrode via a signal generator. 



ions motion in the trap and of the stability region. We 
find Poincarc's sections corresponding to different stabi- 
lity areas and relate them to the observed beam dyna- 
mics. 

The design and operation of the EIBT has been des- 
cribed previously in details [l| and a schematic drawing 
of the ion trap is shown in Fig. [TJ The trap consists 
in a set of coaxial cylindrical electrodes roughly equi- 
valent to two spherical mirrors, the electrostatic ana- 
log of a Fabry-Perot interferometer. The configuration 
of the trap is defined by the potentials of five of these 
electrodes {Vi, V2, 1^3, V4, 14}, the others being grounded. 
The length of the trap is 422mm and the inner radius of 
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the electrodes varies from 8mm to 13mm. An oscillating 
potential Vex siD.{u!ext) can be applied to one of the groun- 
ded electrodes, where both Vex (a few volts) and uJex (a 
few MHz) are adjustable. 

The large number of parameters implied in the tuning 
of the EIBT makes changing setup difficult as it requires 
a lengthy trial and error procedure. In a previous ar- 
ticle [10| , we presented a method to obtain an analytical 
formula of the potential inside the EIBT. Using this me- 
thod, we calculate the potential on the axis V{z) as a 
function of the five potentials applied to the electrodes. 
The general solution of the Laplace equation in cylindri- 
cal symmetry, is a series of the form 



V{r,z)^Y. 
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If we limit ourselves to quadratic terms (we shall come 
back on this approximation below), the trajectory of an 
ion in the trap, without any external excitation, is des- 
cribed by the following set of equations : 
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where, in order to have dimensionless equations, we made 
the following substitutions : z 



L ' 



t 




and Vi ^ where L is the half-length of the trap, 

E is the acceleration voltage and t = — . 

With ;z(0) = 0, r(t) 0, §(i = 0) = 1, Eq. ^ des- 
cribes the motion z{t) along the z-axis, which is periodic 
of period T = l/fz- Substituting in Eq. ((2b)) . we obtain 
the following Hill's equation : 
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This corresponds to a change from the trap reference 
frame to the ion one. The longitudinal motion of the ion 
plays here the same role as the quadrupole traps radio- 
frequency. The principal matrix of (|3]) is 



M{t) 
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where ipi (i, to) is the solution of Eq. ^ with initial condi- 
tions Tpiitoito) ~ 1 and ■0i(to,io) = 0, and '!/'2(t,io) 
with ip2(to,tQ) = and 'ip2ito,to) = 1 respectively. Liou- 
villc's formula [l^l shows that dot M(t, to) — 1 and there- 
fore the characteristic equation of the monodromy matrix 
M {to + T) is given by - 2Aa; + 1 = where : 
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Applying Floquet's theorem [1^, we know that if > 1, 
one of the two solutions is unbound, while for < 1 
there are two bounded solutions : 
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where p±{t + T) = P±{t) and A = cos(7r/3). /3 is often 
called the stability parameter Hill's equation is stable, 
and thus trapping can be observed, when |A| < 1, which 
is equivalent to + 2fc7r < ^ < 1 + 2fc7r Vfc £ Z. 

We performed experiments using ion beams produced 
by the SIMPA[i3| 14.5GHz electron cyclotron resonance 
ion source (ECRIS), accelerated to i? = 5.2keV/charge. 
A magnetic dipole enables us to select the ions by their 
mass/charge ratio. An ion beam is continuously injec- 
ted into the trap, in a potential configuration correspon- 
ding to no entrance mirror and a closed output mirror. 
A sequencer is used to activate a set of fast high- voltage 
switches to close the entrance mirror, trapping the ions, 
which then oscillate between the two mirrors with a per- 
iod of « 2/xs.The number of ions remaining in the trap 
after a fixed time (typically 500/xs) can be measured by 
switching off the exit mirror : particles are ejected out of 
the trap and hit a faraday cup. The signal is processed by 
a charge pre-amplifier and a double delay line amplifier. 
After calibration, the uncertainty on the number of ions 
is ~10%. In this article, we present experimental data 
that was observed with iVe^+ at 5.2kV/charge. We va- 
ried only two of the five potentials, Vi and T4, the three 
others, V2, Vis and V4 being fixed at a constant value, 
respectively 5.85kV, 4.15kV and 1.65kV. 

Figure [5] shows a comparison between the theory out- 
lined above and experiment. The contours show constant 
values |A| — 1 and the stability region is defined by 
A| — 1 < 0. The shaded map represents the number of 
ions remaining after 500/xs of trapping in a given configu- 
ration of (14,14)- For each point, we made several trap- 
ping cycles to obtain the average number of ions trapped 
in this configuration. The resolution of the experimental 
map is lOV along Vz and 50V along Vi. 

Configurations where trapping can be observed expe- 
rimentally are contained in the stability region defined 
by Floquet's theory. However, as we can see on Fig. [51 
some settings (e.g., in the region marked C) are predicted 
theoretically to be stable, which is not observed experi- 
mentally. Complementary to our analytical approach, we 
have developed a 3D C++ code based on the GSL Q 
library, able to simulate individual particle motions in 
the trap, with accurate conservation of the energy and 
momentum. Since the potential is derived from an analy- 
tical formula, these simulations are two order of magni- 
tude faster than what can be achieved with commercial 
software and energy conservation can be guaranteed over 
trapping times of milliseconds. 

The Poincare sections corresponding to the three 
points denoted A, B and C on Fig. [21 evaluated with our 
code, are plotted on Fig. [31 We have taken into account 
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Figure 2. (Color online) Comparison between theory and ex- 
periment : contour lines represent constant values of the stabi- 
lity coefficient | A| — 1 and the shaded map represents number 
of ions remaining after 500/is of trapping in this configura- 
tion. Signals recorded from the charge amplifier connected to 
the faraday cup are presented in the lower part of the figure, 
together with the color code used in the upper part. 




Figure 3. Three Poincare sections corresponding respectively 
to point A, B and C on Fig. [5] The abscissa represents the 
radial position while the ordinate is the radial momentum of 
the particles (dimentionless). We observe a transition from a 
mode centered on r = to a mode centered on r 7^ 0. 



the quadratic radial term in Eq. (|2ap . which leads to a 
system of coupled non-linear equations. We observe that 
point B is located at the transition between a case where 
trapping is concentrated at the center of the trap and a 
case where trapping is centered on a given radius, which 
grows as Vi increases and where chaotic motion occurs. 
This proves that trapping is possible in all the predic- 
ted areas even if the chaotic region is more difficult to 
observe. 



The results above where obtained with no RF exci- 
tation {Vex = 0). We now study the EIBT behavior 
when varying the frequency of the sinusoidal excitation 
signal for Vex = 101^ with a fixed set of electrode po- 
tentials (Vi =7.6kV, V:, =3.23kV). When the excitation 
frequency is resonant with the ions motion frequency, ions 
gain energy at each oscillation and their trajectory be- 
comes unstable as shown in [l5| . These resonances occur 
when 



nfz + mfr = fex, V(r7,, m, k) e I? 



(6) 



where fex — Wex/Stt, = -^Pfz is the modulation of the 
radial motion [see Eq. dS])]. Figure |3] shows that the reso- 
nances appears at the frequencies predicted theoretically 
by the calculation of /3, providing another confirmation 
of this approach. Note that only one species (A^e^^) is 
trapped. This kind of resonance can not only be used to 
monitor the number of ions in the trap, but also leads to 
a new mass spectrometry method. 

The EIBT has already been used as TOF mass spec- 
trometer Q . This technique can be compared to what is 
achieved with Penning traps, where ions oscillations are 
recorded to be analyzed by Fourier transform, giving a 
mass spectrum [l^ . This is only possible if the potentials 
are chosen in order to be in the "synchronization-mode" 
[13, where the ions stay bunched, resisting the coulomb 
repulsion. This counterintuitive phenomenon is explai- 
ned in and requires a minimum number of trapped 
ions (few millions). 

Synchronization restricts the use of the EIBT as a 
mass spectrometer in two ways. First, the unavoidable 
ion losses during trapping limit the time of synchroniza- 
tion and thus of observation, leading to a reduced up- 
per bound of the mass spectrum's resolution. Second, as 
mentioned in Ref. Q, synchronization tends to aggre- 
gate in the same bunch species whose charge/mass ratio 
are close. This is known in mass spectrometry as "peak- 
coalescence" and limits the resolving power. 

The results presented above suggest that an EIBT 
can be used for mass spectrometry in the same way as 
quadrupole traps, where each ion species is successi vely 
ejected from the trap with a parametric excitation [16| . 
Usual methods (e.g., electrospray) could be used to ionize 
samples and send the resulting ions into an EIBT. Since 
the oscillation frequency is proportional to the square 
root of the charge/mass ratio, sweeping the excitation 
frequency fex will provide a mass spectrum of the desired 
range. Besides the simphcity of the setup, this technique 
has two fundamental advantages. First it is independent 
of synchronization and thus is limited neither by synchro- 
nization time nor by peak-coalescence. Second, the stabi- 
lity diagram shown on Fig. [5] is computed only with Eq. 
(121) , whose only parameters are the ratios of the trapping 
potentials to the accelerating potential Vi/ E. This means 
that all the charged particles produced by ionization will 
be trapped independently of their mass and charge. This 
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Figure 4. Experimental observation of the parametric reso- 
nances when fex = nfz+mfr with n and m integers for Ne^'^ 
at 5.2keV/charge. As on Fig. [21 each point represents the ave- 
rage of the ion signal on the faraday cup over several trapping 
cycles of 500/is for each frequency. The peak marked with * 
are not yet understood. 



is very different from quadrupole traps where the stabi- 
lity of the trapping depends on the charge/mass ratio, 
limiting the range of the species that can be trapped and 
analyzed. 

The resonances illustrated on Fig. H] are broad and the 
resolution achieved {M/6M ~ 100) cannot compete with 
advanced techniques. However, many parameters can be 
optimized and all the techniques dedicated to the qua- 
drupole mass filters could be applied to enhance the re- 
solution. For instance, we have noticed that the peaks 
are narrower when Vz is tuned or when the the radial 
eytent nf the injectcH beam is rcHiiccH by closing a, set nf 
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